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Abstract 

This paper deals with the performance analysis of static compensator (STATCOM) based voltage regulator for 
self excited induction generators (SEIGs) supplying balanced/unbalanced and linear/non linear loads. A three- 
phase insulated gate bipolar transistor (IGBT) based current controlled voltage source inverter (CC-VSI) 
known as STATCOM is used for harmonic elimination. It also provides the required reactive power SEIG needs 
to maintain a constant terminal voltage under varying loads. A set of voltage regulators are designed and their 
performance is simulated using SIMULINK to demonstrate their capabilities as a voltage regulator, a harmonic 
eliminator, a load balancer and a neutral current compensator. It also discusses the merits and demerits, to 
select a suitable topology of the voltage regulator according to self excited induction generator. The simulated 
results show that by using a STATCOM based voltage regulator the SEIG terminal voltage can be maintained 
constant and free from harmonics under linear/non linear and balanced/unbalanced loads 

KEYWORDS: Self-excited induction generator, static compensator, voltage regulation, load balancing. 

I. Introduction 

The rapid depletion and the increased cost of conventional fuels have given a thrust to the research on 
self excited induction generator as alternative power sources driven by various prime movers based on 
nonconventional energy sources[5]. These energy conversion systems are based on constant speed 
prime movers, constant power prime movers and variable power prime movers [6] [15]. In constant 
speed prime movers (biogas, biomass, biodiesel etc) based generating systems; the speed of the 
turbine is almost constant therefore the frequency of the generated voltage remains constant. An 
externally driven induction machine operates as a self-excited induction generator (SEIG), with its 
excitation requirements being met by a capacitor bank connected across its terminals. The SEIG has 
advantages [1] [12] [16] [25] like simplicity, being maintenance free, absence of DC, being brushless, 
etc. as compared to a conventional synchronous generator[8][ll][13]. A major disadvantage of an 
SEIG is its poor voltage regulation [14] [24] [18]. It requires a variable capacitance bank to maintain 
constant terminal voltage under varying loads. 

Attempts have been made to maintain constant terminal voltage using fixed capacitor and thyristor 
controlled reactors (TCR), saturable-core reactors and short-shunt connections [6][9][19][21]. The 
voltage regulation provided by these schemes is discrete but these inject harmonics into the generating 
system. However, with the invention of solid state commutating devices, it is possible to make a 
static, noiseless voltage regulator which is able to regulate continuously variable reactive power to 
keep the terminal voltage of an SEIG constant under varying loads. This system, called STATCOM, 
has specific benefits compared to conventional SVC 's [2] [23] [17]. 

Basic topology of STATCOM consists of a 3-phase current controlled voltage source converter (VSC) 
and an electrolytic capacitor at its DC bus. The DC bus capacitor is used to self support a DC bus of 
STATCOM and takes very small active power from SEIG for its internal losses to provide sufficient 
reactive power as per requirements [3] [10]. Here STATCOM is a source of leading or lagging current 
and can be designed in such a way to maintain constant voltage across the SEIG terminals with 
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varying loads. In this paper various STATCOM based VR topologies are presented which are based 
on two leg VSC, three leg VSC for three phase three wire SEIG system[4][7][20]. 
An SEIG is an isolated system, which is small in size, and the injected harmonics pollute the 
generated voltage. The STATCOM eliminates the harmonics, provides load balancing and supplies 
the required reactive power to the load and the generator. In this paper, the authors present a simple 
mathematical model for the transient analysis of the SEIG-STATCOM system under 
balanced/unbalanced. Simulated results show that the SEIG-STATCOM system behaves as an ideal 
generating system under these conditions. 

The brief description about this paper includes, Section2 discusses mainly about the various 
STATCOM controllers used in this paper with the diagrams. Section 3 includes the design of various 
STATCOM techniques included in this paper with the controlling strategies. Section 4 discusses the 
results obtained from the MATLAB/SIMULINK models for various STATCOM techniques applied 
to a self excited induction generator connected to a grid.Section 5 gives the conclusions of this paper. 
The system we tested has the following components: 

• a wind turbine 

• a three-phase, 3-hp,slip ring induction generator driven by the wind turbine 

• various sets of capacitors at stator terminals to provide reactive power to the induction 
generator 

• a three-phase various STATCOM devices 

• a three phase balanced/unbalanced grid 

II. System STATCOM Controllers 

The VRs are classified as three phase three wire VRs and three phase four wire VRs. These VRs are 

based on the two leg VSC, three leg VSC, four leg VSC, three single phase VSC, three leg with 

midpoint capacitor based VSC and transformer based VRs. In the following section, detailed system 

description is presented for different STATCOM based voltage regulators. 

2.1. Three Phase 3- wire voltage regulators 

Mainly two types of VR topologies are discussed for three phase 3 -wire self excited induction 

generator (SEIG). The first one is based on three leg voltage source converter (VSC) and another one 

is based on a two leg VSC with midpoint capacitor. 

2.1.1. Two Leg Voltage Source Converter (VSC) Based Voltage Regulator 

Figure 1 shows an isolated generating system which consists of a constant speed wind turbine, and 

self excited induction generator along with two leg VSC based VR Two legs of VSC are connected to 

each phase of the generator through interfacing inductors while the third phase of the generating 

system is connected to the midpoint of the capacitors. Midpoint capacitors require equal voltage 

distribution across both the capacitors and voltage rating at the DC link of the VSC is comparatively 

higher than the 3 -leg VSC based topology. However switch counts are reduced in this topology of VR 




Figl: Two leg VSC based VR for SEIG system feeding three phase three wire loads. 

2.1.2 Three Leg Voltage Source Converter (VSC) Based Voltage Regulator 

Figure 2 shows an asynchronous generator system based isolated generating system along with three 
leg VSC based STATCOM based voltage regulator. The VR consists of a three-leg IGBT (Insulated 
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Gate Bipolar Junction Transistor) based current controlled voltage source converter, DC bus capacitor 
and AC inductors. The output of the VSC is connected through the AC filtering inductor to the SEIG 
terminals. The DC bus capacitor is used as an energy storage device and provides self-supporting DC 
bus of VR This DC side capacitor supplies the real power difference between the load and SEIG 
during the transient period. In the steady state the real power supplied by the SEIG should be equal to 
the real power demand of the load plus a small power to compensate for the losses of the VR.Thus 
DC capacitor voltage is maintained at a reference value for its satisfactory operation. 
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Fig2: Three leg VSC based VR for SEIG system feeding three phase three wire loads. 

III. Modeling Of SEIG-STATCOM System 

The mathematical model of the SEIG-STATCOM system contains the modelling of an SEIG and 

STATCOM as follows. 

3.1. Modeling of control scheme of STATCOM 

Different components of the SEIG-STATCOM system shown in Fig. are modelled as follows. 

3.1.1. Control scheme of Two Leg Voltage Source Converter (VSC) Based Voltage 
Regulator 

The block diagram of control scheme for two leg VSC based voltage regulator for a SEIG system is as 
shown in Fig. 3 

The control strategy of the two-leg voltage controller based VR is realized similar to three-leg VSC, 
through derivation of reference source currents (i sa r , i sb r ) while main difference between two topology 
to derivation of active component of current as shown in Figure 3. Reference source currents consist 
of two components one is in phase or active power component (i da r , i db r ) for the self supporting DC bus 
of VSC while the other one is in quadrature or reactive power component (i qa r , i qb r ,) for regulating the 
terminal voltage. The amplitude of active power component of the source current (I dm ) is estimated 
using two PI controllers among which, one is used to control the voltage of DC bus of VSC while 
another one is used for equal voltage distribution across the midpoint DC bus capacitors. The output 
of the first PI controller is estimated by comparing the reference DC bus voltage (V dcre f) with the 
sensed DC bus voltage (V dx ).The output of the second PI controller is estimated by comparing the 
voltages across both capacitors (V) and (Va). This voltage error signal is processed using this second 
PI controller The sum of output of both PI controllers (I dm i) and (I dm 2) gives the active power current 
component (I dm ) of the reference source current. The multiplication of I dm with in phase unit amplitude 
templates (u ad ,u bd ) yields the in-phase component of instantaneous reference source currents. These 
(u ad ,u bd ) templates are sinusoidal functions, which are derived by unit templates of in-phase with line 
voltages (u ab ,u bc ,u ca ). These templates (u ab ,u bc ,u ca ) are derived by dividing the AC voltages 
V ab ,V bc ,V ca by their amplitude V t . To generate the quadrature contponent of reference source currents, 
another set of sinusoidal quadrature unity amplitude templates (u aq „u bq , u cq ) is obtained from in-phase 
unit templates (u abd ,u bcd ,u cad ). The multiplication of these components (u aq ,u bq ) with output of the PI 
(Proportional Integral) AC voltage controller (I qm ) gives the quadrature, or reactive power component 
of reference source currents. The sum of instantaneous quadrature and inphase component of source 
currents is the reference source currents (i sa r 4sb r ) and each phase source current is compared with 
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thecorresponding reference source current to generate the PWM switching signal for VSC of the 
controller. 
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Fig.3. Block diagram of control scheme for Two leg VSC based voltage regulator for a SEIG system 

3.1.1.1. Design of Two Leg Voltage Source Converter (VSC) Based Voltage Regulator 

This section presents the detailed design of two-leg VSC based VRs for a SEIG driven by a constant 
speed wind turbine. The two leg VSC and its voltage waveforms are shown in Figure 4. 
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FigATwo leg VSC 

The design procedure is focused on, to determine the value of interfacing inductors, DC link 
capacitors and the voltage across the DC link capacitors along with the rating of the devices. The 
design of the inductor and capacitor depends upon the voltage and current ripples. 

3.1.1.2. Design of the Interfacing Inductor 

In PWM switching of the converter, V C ontroiA,V contro iB and V con troic can be assumed to be constant to be 
constant during one switching frequency time period.At the zero crossing of V contro iA therefore, 



Vc 



^=0 



V C ont r oiB=m a V tri sin(120)= m a y 



(2) 
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Vcontroic=m a V tri sin(240)=- m a y 

When m a is the modulation index and the converter AC terminal voltage vector is defined from the 
line to neutral voltages V an ,V bn ,V cn which can be calculated as follows: 



V an - V aN -V Nn -V aN --(V aN +V bN+ V cN ) 

1 

Vbn = VbN - V Nn= VbN-~(V aN +VbN+V c N) 



V cn - V cN - V Nn= V cN --(V aN +VbN+V cN ) 



(3) 



Where V aN ,V b N and V cN are the converter pole voltages against the midpoint of the DC capacitor and 
V Nn is the voltage between the neutral point(n) and the midpoint of the DC capacitor(N,C). 



Peak to peak inductor current ripple is 

T s 



Cripple = " J n 2 Kimdt=— 

'-•am. u t- 1 ; 



Area A 



(4) 



The interfacing inductor L an ,L bn and L cn can be calculated as follows: 

V3 



^"=1T(5^7)K? + 1 ) 



(5) 



By substituting the values of all parameters, the value of the inductor can be calculated as given in 
table shown in Fig.5. 
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Fig.5. Calculation and selection of various components of two leg VSC based VR 
3.1.2.1. Design of the midpoint D.C link capacitor and its voltage 

Voltage across each capacitor must be more than the peak voltage for satisfactory PWM control as 

2V2"V 



Vdcl-Vdc2- 



V3n 



(6) 



Where m a is the modulation index normally with maximum value 1 .The current which flows through 
the phase connected to the midpoint capacitor is equal to the flow through the capacitors. Therefore 
the ripple in the capacitor voltage can be estimated as: 



' del -ripple - r 
^d 



= rMicdt = 
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(7) 
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^avg ~ \J.y Is 



(8) 



Where I avg is the average current which flows through the DC bus capacitor (C dc i) and I s is the required 
rms compensator current rating of the devices. The voltage rating (V sw ) of the device can be 
calculated under dynamic condition as: 



V sw = (V dc +V d ) 
Where V d is 10% overshoot in the DC link voltage under dynamic conditions. 



(9) 



Rated current which flows through the two leg VSC is I s .The peak value of the current flows through 
the VR considering the safety factor of 1.25 the maximum device current can be calculated as: 



L 



1.25(i 



ripple(pp)" r "-'-s(peak)7 



(10) 



From this voltage(V sw ) and current rating(I sw ) of the IGBT switches can be estimated. Here one design 
example for a two-leg VSC based VR is carried out for feeding 0.8p.f lagging reactive load,SEIG 
requires reactive power of 140-160% of rated generated power. Therefore the additional reactive 
power required from no load to full load at 0.8 lagging p.f load and it is calculated as: 



Additional VAR(Q AR )=V3VI S 
Where V is SEIG line voltage and I s is VR line current. 



(11) 



3.1.2. Control scheme of Three Leg Voltage Source Converter (VSC) Based Voltage 
Regulator 

The block diagram of control scheme for two leg VSC based voltage regulator for a SEIG system is as 
shown in Fig. 6. 
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Fig.6. Block diagram of control scheme for three leg VSC based voltage regulator for a SEIG system 

3.1.2.1. Modelling of control scheme of STATCOM 

Different components of the SEIG-STATCOM system shown in Fig.6 are modelled as follows. 
From the three-phase voltages at the SEIG terminals (V a ,V b and V c ), their amplitude (V t ) is computed 
as: 



The unit vector in phase with V a ,V b and V c are derived as: 

v b 



V a V b 



Uc=" 



v t 



(12) 
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The unit vectors in quadrature with V a ,V b and V c may be derived using a quadrature transformation of 
the in-phase unit vectors u a ,u b and u c as: 



_u b u c 



V3 (u b -u c ) 

W b = — U a H ;=— 

D 2 a 2^3 



W c = 



V3 



U a + 



(Ub-U C ) 

2^3 



"> 



> 



(14) 



3.1.2.2. Quadrature component of reference source currents 

The AC voltage error V er ( n ) at the n th sampling instant is: 



V er (n) - V tre f( n ) - V, 



't(n) 



(15) 



Where V tre f( n ) is the amplitude of the reference AC terminal voltage and V t(n ) is the amplitude of the 
sensed three-phase AC voltage at the SEIG terminals at the n th instant. The output of the PI controller 
(I smq(n)) is used for maintaining constant AC terminal voltage at the n th sampling instant. 
The qudrature components of the reference source currents are computed as: 



saq 



— 1 smq w a \ 1 sbq — 1 smq w b '■> l scq — 1 smq w . 



(16) 



3.1.2.3.In-phase component of reference source currents 



The error in the DC bus voltage of the STATCOM (V dcer(n) ) at the nth sampling instant is: 

Vdcr(n) = Vdcref(n) ~ Vd c (n) (IV) 

where V dcre f(n) is the reference DC voltage and V dc(n ) is the sensed DC link voltage of the STATCOM. 
The output of the PI controller is used for maintaining the DC bus voltage of the STATCOM at the n th 
sampling instant. 
The in-phase components of the reference source currents are computed as: 



1 sad — A smd^a 
1 sbd = I smdU b 
1 sad — -^ smd^c 



(18) 



3.1.2.4. Total reference source currents 

The total reference source currents are the sum of the in phase and quadrature components of the 

reference source currents as :. 

.* .* .* *\ 

+ 1 sad 



1 sa — 1 saq 

1 sb = 1 sbq + 1 sbd 



1 sc — 1 scq "■" 1 scd 



> 



(19) 



3.1.2.5. PWM current controller 



The total reference currents (i sa , i sb and i sc ) are compared with the sensed source currents (i sa , i sb and 
i sc ).The ON/OFF switching patterns of the gate drive signals to the IGBTs are generated from the 
PWM current controller. The current errors are computed as: 



Isaerr — 1 



'Is; 



Isberr — 1 sb Isb 
Iscerr — 1 sc lsc 



(20) 
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These current error signals are amplified and then compared with the triangular carrier wave. If the 
amplified current error signal is greater than the triangular wave signal switch S 4 (lower device) is ON 
and switch Si (upper device) is OFF, and the value of the switching function SA is set to 0. If the 
amplified current error signal corresponding to i saer r is less than the triangular wave signal, switch Si is 
ON and switch S 4 is OFF, and the value of SA is set to 1. Similar logic applies to the other phases. 

3.1.2.6.Modeling of STATCOM 

The STATCOM is a current controlled VSI and is modeled as follows: 

The derivative of its DC bus voltage is defined as: 

i ca SA+i cb SB+i cc SC 

pVdc = (21) 

Where SA, SB and SC are the switching functions for the ON/OFF positions of the VSI bridge 
switches S r S 6 .The DC bus voltage reflects the output of the inverter in the form of the three-phase 
PWM AC line voltage e ab , e bc and e ca . These voltages may be expressed as: 



e ab = V dc (SA-SB) 
e bc = V dc (SB-SC) 

e ca = V dc (SC-SA) 



(22) 



The volt-amp equations for the output of the voltage source inverter (STATCOM) are: 

V a = Rfi ca + Lf pica + e ab -Rfi C b -Lfpicb 1 

V b = R f i cb + L f pi cb + e bc -R f i cc -L f pi cc J (23) 

i ca + icb + ice =0 (24) 

The value of ice from above eqn (24) is substituted into eqn.(23) which results in: 
V b = Rficb + L f pi cb + e bc +R f i ca +L f pi ca +R f i cb +Lfpi cb (25) 

Rearranging the equation results in: 

Lfpica - Lfpicb = V a -e ab -R f i ca +R f i cb (26) 

L f pi ca +2L f pi cb = V b -e bc -Rfica-2R f i cb (27) 

Hence, the STATCOM current derivatives are obtainedby solving eqns. and as: 

• {(Vb-ebc) + 2(Va-e a b)-3Rfi ca } r0Q , 

P lca = ^ (28) 

• {(Vb-ebc)-(Va-e a b)-3Rfi C b} / 00 \ 
P lcb = ^ (29) 

3.1.2.7. AC Line Voltage at the Point of CommonCoupling 

Direct and quadrature axis currents of the SEIG (i ds and i qs ) are converted into three-phases (a, b and 
c). The derivative of the AC terminal voltage of the SEIG is defined as: 



~\T _ {( i a-ira- i ca)-( i b- i rb _i cb)} 



3C 



> (30) 



_ A 7 {(ia-Wc-ica) + 2 (ib-Wb-icb)} 

pv b = 



3C 
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V a +V b +V c = (31) 

Where i a ,ib and i c are the SEIG stator line currents, i ra ,irb and i rc are the three phase load line currents 
and i ca ,icb an d ice are the STATCOM currents. C is the per phase excitation capacitor, which is 
connected across the SEIG terminals. 



IV. Results and Discussions 

The SEIG-STATCOM system feeding linear/non-linear and balanced/unbalanced loads are simulated 
and results are shown in Figs.7-8 . For this study, a 3.5 kW, 440VJ.5A, 4-pole machine was used as a 
generator and the parameters of the generator are given in the Appendix. 

4.1. Performance of two Leg voltage Regulator for a SEIG System 

Here performance of two leg voltage source converter with mid point capacitor based VR topology 
has been simulated using MATLAB/SIMULINK and verified for self excited induction generator 
driven by wind turbine. 
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Fig. 7. Performance of two leg VSC based VR for a SEIG system feeding 3-phase balanced/unbalanced grid 
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Figure 7 shows the transient waveforms of three-phase generator voltages (v abc ), generator currents 

(i abc ), speed, electromagnetic torque, D.C voltage, terminal voltage, load voltage and load current 

respectively demonstrating the response regulating the SEIG terminal voltage while supplying a grid. 

At 0.3 seconds, three-phase nonlinear load is applied and it is found that with application of the sudden 

load, there is increased generator currents, load currents, STATCOM currents and decrease in supply 

voltage due to supplying active and reactive power to the load.The voltage is 75volts and current is 

10A. 

Along with this, short circuit occurs at 0.4seconds,with this there is further increased generator 

currents, load currents, STATCOM currents and decrease in supply voltage due to supplying active 

and reactive power to the load.The voltage is 20volts and current is 15 A. 

At 0.5seconds the STATCOM is connected to the system. Due to this the voltage is reached to the 

required voltage. It is observed that the generator voltage remains constant underbalanced and even 

unbalanced lagging pf loads. Variations in generator speed are observed with the change in load due 

to the drooping characteristic of the wind turbine. 

The STATCOM is disconnected from the system at 0.55 seconds after it reaches the required voltage. 

At 0.6 seconds short circuit is removed and at 0.7 seconds is removed from the load. Now the 

machine is working under steady state conditions. 

The total harmonic distortion (THD) of the generator voltage and current for the three-phase balanced 

case are observed.lt is observed that the THD is less than 5%. 

With the application of the three phase nonlinear loads and short circuits it is found that Voltage 

regulator responds in a desirable manner and maintains constant voltage at the generator terminal. 

Along with this, the DC link voltage and voltage across both midpoint capacitors of voltage regulators 

also remain equal and constant. 

The STATCOM eliminates harmonics so that the generator voltages and currents are free from 

harmonics a scan be observed 

4.2. Performance of three Leg voltage Regulator for a SEIG System 

Here performance of three leg voltage source converter with a capacitor based VR topology has been 

simulated and verified for self excited induction generator driven by wind turbine. 

Figure 8 shows the transient waveforms of three-phase generator voltages (v abc ), generator currents 

(i abc ), speed, electromagnetic torque, D.C voltage, terminal voltage, load voltage and load current 

respectively demonstrating the response regulating the SEIG terminal voltage while supplying a grid. 

At 0.3seconds, three-phase nonlinear load is applied and it is found that with application of the sudden 

load, there is increased generator currents, load currents, STATCOM currents and decrease in supply 

voltage due to supplying active and reactive power to the load.The voltage is 75volts and current is 

10A. 

Along with this, short circuit occurs at 0.4seconds, with this there is further increased generator 

currents, load currents, STATCOM currents and decrease in supply voltage due to supplying active 

and reactive power to the load.The voltage is 20volts and current is 15 A. 

At 0.5seconds the STATCOM is connected to the system. Due to this the voltage is reached to the 

required voltage. It is observed that the generator voltage remains constant underbalanced and even 

unbalanced lagging pf loads. Variationsin generator speed are observed with the change in loaddue to 

the drooping characteristic of the wind turbine. 

The D.C voltage obtained here is having fewer ripples compared with the two leg voltage regulator 

for a SEIG system. Due to these transients there is even change in the variation of the speed and the 

torque. 
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Fig. 8. Performance of three leg VSC based VR for a SEIG system feeding 3-phase balanced/unbalanced grid 

V. Conclusions 

A set of VRs have been designed and their performance have been studied for SEIG system. For 
three-phase three-wire SEIG system two topologies of VR have been demonstrated one is based on 
three leg VSC while another one is based on the two leg VSC. A topology which is based on the two 
leg VSC, requires higher voltage rating of the switches and equal voltage distributed DC link, 
however less number of switching devices are required compared to three leg VSC based topology of 
VR. In three phase three wire SEIG system there are a number of configurations of the VRs for a 
three phase four wire SEIG system . It is observed that the developed dynamic model of the three- 
phase SEIG-STATCOM is capable of simulating its performance while feeding linear/non-linear, 
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balanced /unbalanced loads under transient conditions. From these results, it is found that the SEIG 
terminal voltage remains constant and sinusoidal under a three-phase or a single phase rectifier load. 
When a single-phase rectifier load is connected, the STATCOM balances these unbalanced load 
currents so that the generator currents and voltages remain sinusoidal, balanced and constant and, 
thus, STATCOM acts as a load balancer. A rectifier based non-linear load generates harmonics, which 
are also eliminated by STATCOM. Therefore, it is concluded that the STATCOM acts as a voltage 
regulator, a load balancer and a harmonic eliminator. Although different aspects of uncontrolled 
rectifiers have been modelled as non-linear loads here, the developed model can easily be modified to 
simulate a compensating controlled rectifier as a nonlinear load. For future work they may develop 
the various STATCOM techniques by considering the neutral line and can develop the 3-leg and 4-leg 
wire systems. 
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APPENDIX 

1. STATCOM Control Parameters 

L f = 1.2 mH, R f = 0.045 Q and C dc = 4000uF. 
AC voltage PI controller: K pa =0.05, K ia = 0.04. 
DC bus voltage PI controller K pd = 0.7, K id =0.1 
Carrier frequency = 20 kHz 

2. Parameters of Rectifier Load 

Three-phase rectifier L sL =0.1mH, R SL =1H, Rrl = 22Q, and C RL = 470|iF. 
Single-phase rectifier L SL =0.1mH, R SL = 1 H, R RL =75Q and C RL =150Mf 

3. Machine Parameters 

The parameters of the 3.5 kW,440V, 7.5A, 50 Hz,4-pole induction machine are given below. 
R s = 0.69 Q, R r = 0.74Q, L ls = L lr = 1.1 mH, J = 0.23kg/m2, 
Lss = Lis + L m and L rr = L lr + L m 

4. Terminal capacitor 
C = 15 uF/ phase 

5. Air gap voltage: 

The piecewise linearization of magnetization characteristic of machine is given by: 

E 1= " X m >260 

E!=1632.58-6.2X m 233.2<X m <260 

E!=1314.98-4.8X m 214.6<X m <233.2 

E 1 =1183.11-4.22X m 206<X m <214.6 

E!=l 120.4-3.9.2X m 203.5<X m <206 

E!=557.65- 1 . 144X m 197.3<X m <203.5 

E!=320.56-0.578X m X m <197.3 
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